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ABSTRACT 

We announce the discovery of a highly inflated transiting hot Jupiter discovered by the KELT- 
North survey. A global analysis including constraints from isochrones indicates that the V = 10.8 
host star (HD 343246) is a mildly evolved, G dwarf with = STSd+gg K, log 5 = 4.078^Q;g54, 

[Fe/H] = 0.272 ± 0.038, an inferred mass M* = 1.21lto;o66 radius = 1.67to;i2 R-o- The 

planetary companion has mass Mp = 0.867tQ;gg4 Mj, radius Rp = 1.86to;i6 Rj, surface gravity 
\oggp = 2.793lQ;g75, density pp = 0.167to;g38 g cm“^. The planet is on a roughly circular orbit 
with semimajor axis a = 0.04571^Q;gQgg4 AU and eccentricity e = 0.035^Q;g25- The best-fit linear 
ephemeris is Tg = 2456883.4803 ± 0.0007 BJDtdb and P = 3.24406 ± 0.00016 days. This planet is 
one of the most inflated of all known transiting exoplanets, making it one of the few members of a 
class of extremely low density, highly-irradiated gas giants. The low stellar logi/ and large implied 
radius are supported by stellar density constraints from follow-up light curves, plus an evolutionary 
and space motion analysis. We also develop a new technique to extract high precision radial velocities 
from noisy spectra that reduces the observing time needed to confirm transiting planet candidates. 

This planet boasts deep transits of a bright star, a large inferred atmospheric scale height, and a high 
equilibrium temperature of Teq = K, assuming zero albedo and perfect heat redistribution, 

making it one of the best targets for future atmospheric characterization studies. 

Subject headings: planetary systems — stars: individual (HD 343246) 
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1. INTRODUCTION 

Ground-based surveys for transiting exoplanets have 
been extraordinarily successful in detecting Jupiter-size 
planets orbiting very close to their host stars (hot 
Jupiters). These planets are some of the easiest to detect 
and characterize but also some of the most intrinsically 


(IGould et al. 12006 

Gumming et al. 12008 Bayliss & Sack- 

ett 

2011 Wright et al.||2012 Howard et al. 

2012). Both 


radial velocity (KV j surveys and transit surveys are heav¬ 
ily biased to detect these massive close-in objects yet we 
still know of only 162 exampleof hot Jupiters out of 
the thousands of currently known exoplanets. 

Due to their short orbital periods and large transit 
and/or radial velocity (RV) signals, these planets make 
excellent laborat ories to study planet formation and mi- 


Mordasini et al.IzOOy 

), atmospheric properties and com- 

position (jZhao et al. 

2014 

Lockwood et al. 

2014 

Sing 


planets (Snellen et al. 2014). Several mysteries remain 
unsolved tor this population of exoplanets including the 
differences between the orbital distance distribution of 
close-in Jupiters orbiting metal rich and metal poor stars 


(Dawson & Murray-Clay [2013 ), the freq uency of long- 


period companions in hot Jupiter systems (Knutson et al. 


2014| |Ngo et al.|2015|) , the source of large misal ignments 
between the orbital axis and stellar spin axis (Albrecht 
et al.|2012 Dawson]2014 ), and the reason that many hot 
Jupiters "are inflated to extremely large radii (Batygin 
et al.||20111 |Ginzburg & Sari 2015). Each hot Jupiter 
-fiances ou 


discovery enhances our ability to explore these phenom¬ 
ena and study individual systems in exquisite detail. 

RV surveys have discovered many non-transiting hot 
Jupiters around very bright stars, while the target stars 
of transit surveys tend to be much fainter. This is pri¬ 
marily due to the fact that only a narrow range of or¬ 
bital parameters cause a planet to transit and many tens 
of thousands of stars must be observed in order to de¬ 
tect a single transiting planet. Since bright stars are 
distributed across a large area of the sky it is gener¬ 
ally more efficient to observe small areas of the sky to 
a greater depth for transits of fainter stars. The Kilode- 
gree Extrem ely Little Telescope survey (KELT, Pepper 


et al. 


2007) serves to bridge the gap in host star bright- 


ness between RV-detected and transit-detected planets 
by using extremely small aperture telescopes with large 
fields of view to observe nearly the entire sky. The survey 
is optimized to detect planets orbiting stars with V-band 
magnitudes between 8 and 10 but it has also been suc¬ 
cessful i n detecting planets around slightly fainter stars 
as well ( Collins et al.|[2M4 ). 


We report the discovery of a highly inflated (i?p = 
l-86lo;^g — 0-867lQ;Qg® Mj) hot Jupiter orbit¬ 

ing the moderately bright V=10.84 G dwarf; HD 343246 
(KELT-8 hereafter). KELT-8b has the 2nd largest ra¬ 
dius and 7th lowest density among all tra nsiting exo- 


et al.|2012 

Bento et al. 

2014 

). The planet lies well above 

the 

Seager et al. 

12007 

) mass-radius relation tor a pure 


hydrogen composition. In order to facilitate rapid and ef¬ 
ficient RV confirmation of the planet we developed a new 
technique that saves significant telescope resources when 
collecting RV meas urements using the iodine technique 
(|Butler et al.|1996D. In p|we describe our discovery and 
follow-up observations. \Ve describe our spectroscopic 
analysis of the host star and summarize the inferred stel¬ 
lar properties in ^ A close companion stellar object 
that was ultimately foun d to be a background contami¬ 
nant is described in §3.4[ We outline our new technique 
for RV extraction in A global analysis of the photo¬ 
metric and RV data is presented in il and we conclude 
with a discussion in m and summary in 0 

2. OBSERVATIONS 
2.1. KELT Photometry 

The KELT-North survey telescope consists of an 
Apogee AP16E imager (4K x 4K GGD with 9^m pixels) 
and an 80mm Mamiya 645 camera lens (42 mm diameter, 
f/1.9) behind a Kodak Wratten #8 red-pass filter. This 
setup achieves a 26° x 26° field of view with roughly 23” 
pixel“^. For a mo re complete descript ion of hardware 
and operations, see|Pepper et al.| (|2007|). 

Raw KELT science images are dark-subtracted and 
flat-fielded using standard methods and then reduced us¬ 
ing a heavily-modifie d implementation of the ISIS image 
subtraction package ( Alard fc Lupton|1998 Alard|[2000| 
coupled with our own high-performance background- 
removal routines. Stars are identified for extraction using 
the standalone DAOPHOT II PSF-fitting software pack¬ 
age (|Stetson||1987l |I990|). To reduce systematic errors, 
extracted light curve s are processed with the trend filter¬ 


ing algorithm (TFA; Kovacs et al.j|2005 ) prior to period 
search. A more complete description of the data reduc- 
ti on pipeline and cand idate selection process are available 


Siverd et al. (2012). 


lELi’-8 resides in KELT-North field 11, centered on 
{a = 19'*26’"48", 5 = -b31°39'56"; J2000). The data set 
consists of 5978 images acquired between 30 May, 2007 
and 14 June, 2013. The full discovery light curve, phased 
to the KELT-8b transit ephemeris, is shown in Figure 

2.2. Follow-up Photometry 

We acquired follow-up time-series photometry of 
KELT-8 to better determine the transit shape and to 
check for a consistent transit depth across the optical fil¬ 
ter bands. To schedule foll ow-up photom etry, we used 
the Tapir software package ([Jensen 


2013). We obtained 


Planets with orbital periods shorter than 10 days and 
M sin 2 > 0.5 Afj based on a 2015 Mar 3 query of exoplanets.org 


full primary transits in multiple bands between July 2014 
and September 2014. Figure shows all the primary 
transit follow-up light curves assembled. A summary of 
the follow-up photometric observations is shown in Ta¬ 
ble We find consistent ratios in the g', r', i', 

V, and Ic filters after compensating for blending with a 
nearby neighbor in some observations (see below), help¬ 
ing to rule out false positives due to blended eclipsing 
binaries. The lower panel of Figure shows all primary 
transit follow-up light curves from Figure [^combined and 
binned in 5 minute intervals. This combined and binned 
light curve is not used for analysis, but rather to show 
the overall statistical power of the follow-up photometry. 

AstroImagej[^ (AIJ) (K. A. Gollins & J. F. Kielkopf 

http://www.astro.louisville.edu/software/astroimagej/ 
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Fig. 1.— Discovery light curve of KELT-8b from the KELT-North telescope. The photometry is phase-folded to the ephemeris listed in 
Table |4] The red line is the same data binned in intervals of width 0.01 in orbital phase (ps 46.7 minutes). 


2015, in preparation) was used to calculate the differ¬ 
ential photometric data from all calibrated image se¬ 
quences. AIJ is a general purpose image processing pack¬ 
age, but is optimized for processing time-series astronom¬ 
ical image sequences. 

2 . 2 . 1 . MORC 

We observed one complete KELT-8b transit on UT 
2014-08-04 from the 0.6 m Moore Observatory RCOS 
(MORC) telescope, operated by the University of 
Louisville. MORC is equipped with an Apogee U16M 
4K X 4K CCD giving a 26’ x 26’ field of view and 0.39 
arcseconds pixel” The transit was observed pseudo- 
simultaneously in two filters by alternating between the 
g' and i' filters from one exposure to the next. The filter 
change is made during CCD read-out to minimize loss of 
time on the the sky. Since KELT-8b has a nearby neigh¬ 
bor, telescope defocusing was minimal to provide better 
separation of the flux from the two stars on the detector. 

2 . 2 . 2 . GCorfini 

We observed one full transit on UT 2014-07-31 from 
Giorgio Corfini’s home observatory (GCorfini) in Lucca, 
Italy. The observations were obtained in V-band using 
a 0.2 m Newtonian telescope and a SBIG STT-6303 ME 
CCD 1536 X 1024 pixel camera, giving a 59’ x 39’ field 
of view. No defocusing was applied, but the target and 
neighbor star were severely blended on the detector due 
to seeing. A photometric aperture was chosen that in¬ 
cluded both KELT-8 and its neighbor. The resulti^ 
light curve was corrected for blending as described in 

2 . 2 . 3 . ZRO 

We observed one full transit in the Vc filter on UT 
2014-07-31 at Roberto Zambelli’s personal observatory 
(ZRO). The telescope was not defocused and both stars 
were cleanly resolved from each other. The observations 
were obtained using a Meade LX 200 ACE GPS 12” tele¬ 
scope, a f/6.3 focal reducer, and an SBIG ST8 XME 1530 
X 1020 GGD with a pixel scale of 0.92 arc seconds pixel”^ 
and a 23’ x 15’ field of view. 

2 . 2 . 4 . WCO 

We observed a full transit in the r' band at Westmin¬ 
ster College Observatory (WCO) in Pennsylvania on UT 


TABLE 1 

Summary of Photometric Observations 


Telescope 

UT 

(2014) 

# 

Obs 

Filter 

Band 

Cycle^ 

(sec) 

RMSt" 

(10”3) 

PNR'^ 

' min ' 


GCorfini 

7-31 

73 

V 

243 

3.0 

6.0 

1.1 

ZRO 

7-31 

110 

V 

216 

3.5 

6.6 

1.9 

MORC 

8-04 

65 

g' 

127 

2.1 

3.0 

2.2 

MORC 

8-04 

62 

i' 

127 

2.3 

3.3 

2.6 

WCO 

8-30 

209 

r’ 

81 

3.8 

4.4 

3.1 

CROW 

9-09 

49 

Ic 

332 

2.3 

5.4 

1.3 


Cycle time in seconds, calculated as the mean of exposure time plus 
dead time during periods of back-to-back exposures. The MORC g' and 
i' exposures were alternating, so cycle time is calculated for exposures 
in both filters combined. 

RMS of residuals from the best fit model in units of 10“® . 

Photometric noise rate in units of 10“^ minute”^, calculated as 
RMS/\/r, where RMS is the scatter in the light curve residuals and 
P is the mean number of cycles (exposure time and dead ti me! per 
minute dur ing periods of back-to-back exposures (adapted from|Fulton| 
|et al.|201l| |. 

Scaling factor applied to the uncertainties to ensure that the best 
fitting model has !}4- 

2014-08-30. The observations were obtained using a Ce- 
lestron 0.35 m C14 telescope with an SBIG STL-6303E 
3072 X 2048 GGD, giving a 24’ x 16’ field of view and 
0.5 arcseconds pixel”^. The data were collected with the 
telescope in focus and the two stars were resolved from 
each other. 

2 . 2 . 5 . CROW 

We observed one full transit in the Ic band at Canela’s 
Robotic Observatory (GROW) in Portugal on UT 2014- 
09-09. The observations were obtained using a 0.3 m 
LX200 telescope with an SBIG ST-8XME 1530 x 1020 
CCD, giving a 28’ x 19’ field of view and 1.11 arcseconds 
pixel” No telescope defocus was applied, but KELT-8 
and the neighbor were blended on the detector due to 
seeing. A photometric aperture was selected to encircle 
the flux from both stars. The resulting light curve was 
corrected for blending as described in 

2 . 3 . High-resolution Spectroscopy 
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TABLE 2 

Radial Velocities of KELT-8 


BJDtdb 
(- 2440000) 


RV 

(ms-l) 


ilrv 

(ms-l) 


BS 

(ms“^) 


ciBS 

(ms-i) 


16884. 

16886. 

16888. 

16890. 

16891. 

16892. 

16894. 

16895. 

16900. 

16903. 

16905. 

16906. 

16970. 


7240398 

7223760 

6894588 

6920118 

6849308 

6862628 

6776168 

6940738 

7203146 

7467735 

6816215 

7226164 

6342444 


-52.700 

16.089 

54.770 

-85.923 

39.139 

120.081 

-30.453 

98.475 

-99.117 

-104.461 

91.985 

-90.377 

55.520 


13.698 

11.487 

11.598 

11.536 

11.894 

10.732 

11.729 

11.258 

11.290 

12.930 

13.023 

12.995 

20.228 


18.194 

60.348 

21.460 

104.118 

-18.528 

15.611 

-18.657 

18.095 

-3.809 

67.604 

6.130 

117.227 

73.214 


14.832 

22.517 

21.508 

56.715 

26.524 

12.552 

35.822 

18.990 

25.208 

41.392 

15.736 

77.109 

34.466 


strumental PSF and used as the stellar template in the 
forward modeling process. However, in this case the star 
is too faint to collect the signal-to-noise needed for re¬ 
liable deconvolution in a reasonable amount of time on 
the APF. Instead we simulate this observation by using 
the SpecMatch software (Petigura et al. 2015 in prep.) 
to construct a synthetic template from thejCoelho (20141 
models and best fit stellar parameters. I'his process is 
described in more detail in Section |4l 



- 2-10 1 2 
Time - (hrs) 


Fig. 2.— Top: Follow-up photometry of KELT-8b primary tran¬ 
sits. The source, date, and filter for each transit is annotated. The 
best-fit models are shown in red. Bottom: All follow-up transit 
light curves combined and binned in 5 minute intervals. This light 
curve is not used in the analysis and is shown simply to demon¬ 
strate the combined photometric precision of the follow-up light- 
curves. 


We collected a total of 13 RV measurements of KELT- 
8 using the Levy high-resolution optical spectrograph 
mounted at the Nasmyth foci of the Autom ated Planet 


Finder Telescope (APF) at Lick observatory (Vogt et al. 

' L] 20141). I'he mea- 

TT noi 1 


2014 Radovan et al. 2014 Burt et al. 


surements were collected between U'i 2014 August 15 
and UT 2014 November 9 and are presented in Table 
We observed this star using the 2 x 8" slit for a spectm 
resolution of R « 80,000. We pass the starlight through 
a cell of gaseous iodine which serves as a simultaneous 
wavelength and point s pread function (PSF) reference 
(Marcy & Butler 1992). Relative radial velocities are 
calculated by tracing the doppler shift of the stellar spec¬ 
trum with respect to the dense forest of iodine lines using 
a forw ard modeling technique described in [Butler et al.| 
(1996). Traditionally, a high signal-to-noise iodine-free 
observation of the same star is deconvolved with the in- 


2.4. Adaptive Optics Imaging 
We acquired visible-light adapti ve optics images of 


KELT -8 using the Robo-AO system (Baranec et al. 12013 


20141 on the 60-inch Telescope at Palomar Observatory. 
(Jn UT 2015 March 8, we observed KELT-8 in the Sloan- 
i' filter as a sequence of full-frame-transfer detector read¬ 
outs at the maximum rate of 8.6 Hz for a total of 120 
s of integration time. The individual images are cor¬ 
rected for detector bias and flat-fielding effects before 
being combined using post-facto shift-and-add process¬ 
ing using KELT-8 as the tip-tilt star with 1 00% frame 
select ion to synthesize a long-exposure image (Law et al. 

Mil. 


We calculate the 5 a contrast limit as a function of 
angular separation by first determining the background 
noise level for concentric rings of width equal to the 
FWHM moving outward from the primary star. A simu¬ 
lated companion (a dimmed cutout of the primary star) 
is inserted into the closest ring, with a random PA. 
An auto-companion detection code then searches for the 
simulated companion. If it is found, the companion is 
dimmed further and reinserted. This is continued un¬ 
til the companion is not found with a confidence greater 
than 5 a. We repeat this process for each annulus, and 
fit the sparse measurements with a function of the form 
a —b/(r — c), where a, b, and c are free parameters in the 
fit and r is the radius from the target star. We convert 
contrast lim its to mass limits using the models of Baraffe 
et al. (2002|). The resulting contrast and mass limits are 
presented m Figure]^ 


3. STELLAR PROPERTIES 
3.1. Keck/HIRES Spectroscopy 

In order to obtain precise values for the stellar param¬ 
eters we collected a moderate signal-to-noise iodine-free 
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Fig. 3.— Top: Robo-AO image of KELT-8 showing the back¬ 
ground giant star 8^8 to the SW of KELT-8. Bottom: Five a 
contrast limit near KELT-8 derived from the Robo-AO image. 


observation using the HIRES s pectrograph on the Keck 
I telescope (Vogt et al. 1994). We measured the ef¬ 
fective temperature ('Jes), surface gravity (log 5 ), iron 
abundance ([Fe/H]), and rotational velocity of the star 
(usini) using the tools available in the SpecMatch soft¬ 
ware package (Petigura et al. 2015 in prep). We first 
correct the observed wavelengths to be in the observer’s 
rest frame by cross-correlating a solar model with the 
observed spectrum. Then we fit for Teg, logg, [Fe/H], 
usini, and the instrumental PSF using the underlying 
different ial-evolution Mar kov Chain Monte Carl o (DE- 
MCMC, Ter Braak 2006) machinery of ExoPy (Fulton 


et al. 20l^). At each point in the MCMC chains a syn- 
thetic spectrum is created by interpolating the [Coelho 


(2014) grid of stellar models for a set of Teff, log g, [Fe/HJ 
values and solar alpha abundance. We convolve this syn¬ 
thetic spectrum with a rotational plus mactroturbulance 
broadening ke rnel using the prescript ions of [Valenti fc| 


Fischer (2005) and Hirano et al. (2011). Finally, we per¬ 
form ahother convolution with a gaussian kernel to ac¬ 
count for the instrumental PSF and the synthetic spec¬ 
trum is compared with the observed spectrum using 
to assess the goodness of fit. The priors are uniform in 
Teff, log 5 , and [Fe/H], but we assign a gaussian prior 
to the instrumental PSF that encompasses the typical 
variability in the PSF width caused by seeing changes 
and guiding errors. Five echelle orders of the spectrum 


are fit separately and the resulting posterior distributions 
are combined before taking the median values for each 
parameter. 

Parameter uncertainties are estimated as the scatter 
in spectroscopic parameters given by Spec Match relative 
to the values in Valenti & Fischer (2005) for 352 of the 
stars in their sample and 76 stars in the [Huber et al.j 
(2012) asteroseismic sample. Systematic trends in Spec- 
IVlatch values as a function of Teg, logg, and [Fe/H] rela¬ 
tive to these benchmark samples are fit for and removed 
in the final quoted parameter values. More details about 
the benchmark comparison and calibrations will be avail¬ 
able in Petigura et al. (2015, in prep). These spectro¬ 
scopic parameters are then used as input priors for the 
global fit discussed in The best fit spectroscopically- 
determined stellar parameters are log 5 = 4.23 ± 0.08, 
[Fe/H] = 0.272 ± 0.038, and Teg = 5754+1^ K. 

The log g derived from spectroscopy is only marginally 
consistent with the log 5 derived from the global fit 
(4.078;HS)- The log g value measured from the transit 
light curves is likely more reliable because log g can be a 
difficult quantity to measure from high-resolution spec¬ 
troscopy alone. We perform a global analysis both with 
and without the spectroscopic log 5 as a prior (see ©• 

We also collected a single exposure with the iodine cell 
in the light path with a slightly higher S/N ratio to es¬ 
tablish a long term baseline for possible followup of this 
target with Keck in the future. We used this observation 
to extract the logT^j^ stellar activity metric. We find 
that logi?yp;=-5.108 which indicates that the star very 
chromospherically quiet. 


3 . 2 . Spectral Energy Distribution 

We estimate the distance and reddening to KELT -8 
by fitting the B-band through W ISE W4-band s pectral 
energy distribution (SED) to the Kurucz (1979) stellar 
atmosphere models. We als o found a GALEX i\ U V mea¬ 


surement in the literature ( Bianchi et al.[[ 2 MT ), however 
this measurement shows a large excess in hux relative to 
the SED compatible with all other measurements. Since 
the PSF of GALEX in the NUV band is 4.9” we suspect 
the GALEX measurement is contaminated by the visual 
companion 878 away and we omit this measurement from 
the fit. The WISE measurements also have a large PSF 
relative to the separation between the two components, 
but we do not find any evidence for significant contam¬ 
ination. We fix Teff, logg, and [Fe/H] to the best fit 
values in Table and leave distance (d) and reddening 
(Ay) as free parameters. We find best-fit parameters of 
Ay = 0.15 ± 0.06 magnitudes and d = 236 ± 9 pc with a 
reduced = 3.7 (considering only the flux uncertainty 
and not the uncertainty in the model). This a grees quite 
well with the d istance of 233 pc estimated by Pickles <£] 
Depagne (2010). 


3 . 3 . Evolutionary State 

We estimate the age of KELT -8 by fitting Yonsei-Yale 
isochrones to the values of Teg, logg, and [Fe/H] given 
in Table [4| We fix the stellar mass to the value of 
M* = 1.211 Mq listed in Table Our best fit stellar 
parameters indicate that the star happens to fall on a 
very rapid part of the evolutionary track, the so-called 
“blue hook” just prior to crossing the Hertzsprung gap 
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Fig. 4. — SED fit to KELT-8. The red crosses are the photom¬ 
etry measurements. The vertical errors are the 1 cr photometric 
uncertainties, and the horizontal error bars are the effective widths 
of the passbands. The leftmost point (GALEX NUV) was omitted 
from the fit due to possible contamination from the neighboring 
star at 8”8 separation. The black line is the best-fit stellar at¬ 
mosphere model and the blue points are the predicted passband- 
integrated fluxes of the model corresponding to our observed pho¬ 
tometric bands. 

(see Figure]^. In principle, this would allow for a very 
precise measurement of the stellar age, however, the mea¬ 
surement errors on T^s and log g don’t allow us to locate 
the star’s exact position within the blue hook. We con¬ 
servatively quote an age range of 4.9-5.8 Gyr that spans 
the entire blue hook. We note that the uncertainty on the 
age likely does not follow a Gaussian distribution because 
the rapid evolution through the blue hook means that 
stars on the blue hook should be rare.It is much more 
likely that it falls on one side or the other of the blue 
hook simply because stars spend such a small fraction 
of the lifetimes on the blue hook. A more detailed evo¬ 
lutionary analysis that includes rotation, chemical com¬ 
position, and a full consideration of all relevant priors 
is needed to determine the exact evolutionary state of 
KELT-8. 



TefftK] 


Fig. 5.— Yonsei-Yale isochrone fit for KELT-8. The solid black 
line is the isochrone for a star with the mass and metallically 
of KELT-8 listed in Table The red cross is the position of 
KELT-8, and the blue points label various ages along the evolu¬ 
tionary track. The green cross is the position of KELT-8 using the 
spectroscopically-determined log g. The “blue hook” is the kink in 
the track where the star’s Tgfr increases while log g decreases for a 
short amount of time between an age of 4.9 and 5.8 Gyr. 


3.4. Unrelated visual companion 

KELT-8 was identified as a vis ual binary in 19 05 by 
French astronomer Abel Pourteau (|Pourteau|1933|). Fig¬ 
ure 1^ shows KELT-8 and the fainter visual companion 
8"8 to the SW clearly resolved. In addition to the discov¬ 
ery astrometry, further astrometric measur ements were 
made in 1951 (iG ellera 1984 1990), 200 0 (Gutri et al. 
20031, and 20011 Hartkopt et al. 2013b). 'These mea- 
surements all show a trend ot decreasing separation and 
increasing position angle, although the scatter in the pre- 
GCD measurements is large (see Figure]^. We also mea¬ 
sured astrometry from our Robo-AO data and find that it 
is consistent with this trend, especially if only the GGD 
measurements are considered. This is consistent with 
KELT-8 moving with the constant proper motions given 
in Table and the companion being stationary in the 
background. 

We also collected an APF spectrum of the secondary 
component and used the SpecMatch pipeline to extract 
stellar parameters. We cross-correlated the best-fitting 
model with the observed spectrum after correcting for 
the Earth’s barycentric motion to determine an absolute 
systemic velocity. SpecMatch returned T^s = 4702 ± 60, 
log 5 = 2.77 ± 0.08, and [Fe/H] = 0.01 ± 0.04 which 
is consistent with the companion being a G or K giant 
with an intrinsic luminosity much greater than that of 
KELT-8. Since this component appears approximately 
one magnitude fainter than KELT-8 in V band it must 
be in the background. In addition, the systemic RV of 
the secondary of 34 ± 1 kms“^ is highly discrepant with 
the velocity of KELT-8 (-31.2 ± 0.3 kms“^). These lines 
of evidence all suggest that the two stars are physically 
unrelated. 


3.5. UVW Space motion 

We evaluate the motion of KELT-8 through the galaxy. 
The absolute radial velocity measured from the APF 
spectra of KELT-8 is —31.2±0.3 km s“^, and the proper 
motions from Tycho-2 are —13.7 ± 1.4 and —29.5 ± 1.3 
milli-arcseconds per year in RA and DEG respectively. 
These values transform to U,V,W space motions of 
20.6 ± 1.7, —32.6 ± 1.2, 0.7 ± 1.5 km s“^ respectively, 
making KELT-8 an un ambiguous member of the thin 
disk stellar population (Bensby et al. 2003). Since it is 
not a halo star this further supports our adoption of the 
stellar log g derived from the transit fit which gives stel¬ 
lar parameters (particularly met allicity) that are m uch 
more common for thin disk stars (Ivezic et al.||2012 ). 


4. SYNTHETIG TEMPLATE RADIAL VELOCITIES 

Iodine-based radial velocities are traditionally derived 
using a forward modeling process that requires three pri¬ 
mary ingredients; a template spectrum of the star with 
the iodine out of the light path, an ultra-high resolution 
spectrum of the iodine absorption for the particu lar cell 
in use, and a description of the instrumental PSF (Butler 


et al. 1996). The stellar template and the iodine spec¬ 


trum are multiplied together then convolved with the 
PSF to match the observed spectrum taken through io¬ 
dine. The stellar template is usually an observation of the 
same star taken at 2-3 times higher signal-to-noise (S/N) 
ratio than the iodine observations and using a narrower 
slit for higher resolution. Since the PSF is constantly 
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TABLE 3 

Adopted Stellar. Properties oe KELT-8 


Parameter (units) 

Value 

Source 

Reference 

Names 

HD 343246 

SIMBAD 




GSC 02109-00049 

SIMBAD 




2MASS J18531332-I-2407385 

SIMBAD 




TYC 2109-49-1 

SIMBAD 




WDS 18533+2407 

WDS iHartkopf et al.l ||2013a|| 

Spectral type 

G2V 

this work 




RA (J2000) 

18 53 13.3216 

Tycho-2 

H0g et al. 1 

2000 


DEC (J2000) 

+24 07 38.603 

Tycho-2 

H 0 g et al. I 

^UUU 


(kms“l ) 

20.6+ 1.7 

this work 



V (kms“^ ) 

-32.6+ 1.2 

this work 



W (kms-i ) 

0.7+ 1.5 

this work 



RA proper motion (mas/yr) 

-13.7+1.4 

Tycho-2 

H0g et al. I 

2000 


DEC proper motion (mas/yr) 

-29.5+ 1.3 

Tycho-2 

H 0 g et al. 1 

liUUU 


Systemic velocity (kms“^ ) 

-31.2 + 0.3 

this work 



NUV (mag) 

16.884 + O.R 

GALEX i 

hanchi et al 

(12011 1 

Bt (mag) 

11.713 + 0.057 

Tycho-2 

H 0 g et al. I 

20DD 


Vt (mag) 

10.925 + 0.048 

Tycho-2 

ti 0 g et al. 

2DDD 


B (mag) 

11.545 + 0.065 

ASCC 

RnarcnenKo 

1 2UU 


V (mag) 

10.833 + 0.054 

ASCC 

Kharchenko 

2UUi 


J (mag) 

9.586 + 0.026 

2MASS 

Uutri et al. 

12UU3 


H (mag) 

9.269 + 0.032 

2MASS 

Cutri et al. 

I20U3 


K (mag) 

9.177 + 0.021 

2MASS 

Cutri et al. 

1 2UU3 


W1 (mag) 

11.664 + 0.05 

WISE 

Uutri et al. 

(2014 


W2 (mag) 

12.302 + 0.05 

WISE 

Cutri et al. 

i2U14 


W3 (mag) 

14.17 + 0.05 

WISE 

Cutri et al. 

I2U14 


VE4 (mag) 

15.193 + 0.339 

WISE 

Uutri et al. 

(2014 


nsin2(kms~^) 

3.7+ 1.5 

this work 



Age (Gyr) 

5 4+0-4 

this work 



Distance (pc) 

236 + 9 

this work 



Av 

0.15 + 0.06 

this work 



log 

-5.108 

this work 




Positive U is in the direction of the Galactic Center. 
Likely contaminated by the neighbor star 878 away 


changing on a slit-fed spectrograph due to changes in 
the slit illumination caused by variable seeing and guid¬ 
ing patterns, we need to deconvolve this stellar template 
to remove distortions in the line profile. This deconvolu¬ 
tion process also requires a measurement of the instru¬ 
mental PSF. By construction, we do not have a simulta¬ 
neous PSF calibration for the template observations as 
we do for the iodine observations. We rely on brack¬ 
eting observations of rapidly rotating B stars through 
iodine, but the PSF is never exactly the same for any 
two observations and this non-simultaneous calibration 
introduces errors. In addition, the deconvolution process 
is inherently unstable because we are essentially trying 
to create a view of the star at infinite spectral resolution 
using data taken at finite spectral resolution. Collecting 
as much signal as possible for the template observation 
helps ensure that this component of the model is not the 
limiting factor on the final RV precision. 

For faint stars (V>12), iodine exposures on 
Keck/HIRFS can take 45-60 min each. Collecting 
a template at 2-3 times higher S/N ratio would take a 
significant fraction of a night. At APF the situation 
is even worse due to its much smaller aperture (2.4 
m vs. 10 m). When collecting reconnaissance RVs 
for a new candidate from a transit survey to look for 
various false positive scenarios only 2-3 RVs may be 
needed before a false alarm scenario can be confirmed 
or eliminated. Collecting a high S/N ratio template can 
take as much time as the entire set of RV observations, 
effectively doubling the observational time needed in 


these scenarios. 

In order to circumvent these problems we developed 
a techn i que t hat uses a synthetic template using the 


Coelho (20141 models and spectral synthesis code from 
SpeciVIatch. We first determine the spectral parameters 
from the best available information. Approximate values 
for Teff, logg, and [Fe/H] can be obtained using spectral 
regions outside of the iodine absorption regime, photom¬ 
etry, or low resolution spectroscopy. We normally use a 
single echelle order redward of the iodine regime to esti¬ 
mate usini. Tgfj and usini have the most drastic impact 
on the resulting RV, but fortunately Teff is the easiest 
parameter to estimate in the absence of a high resolu¬ 
tion spectrum and v sin i can be estimated with very few 
spectral lines. Once the spectral parameters are estab¬ 
lished we synthesize the full spectrum within the iodine 
region setting the instrumental PSF width to zero. This 
provides an extremely high-resolution, noise-free stellar 
template to use as one of the inputs to the normal RV ex¬ 
traction code. This foregoes the need for deconvolution 
and saves valuable telescope time. 

This technique is not perfect due to inaccuracies in 
the model spectra. Spectral lines are missing in some re¬ 
gions, blended lines may only be accounted for as a single 
line, or the central wavelengths of lines may be slightly 
inaccurate. This means that using these synthetic tem¬ 
plates we can never achieve the 2-3 ms“^ precision that 
we can obtain using the traditional technique for bright 
stars. Figure shows velocity RMS for the RV standard 
star (HD 940^ as a function of S/N ratio of the iodine 
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Fig. 6. — Astrometry measurements of KELT-8. The green 
squares were measured using astrometric eyepieces or photographic 
plates. The black points were measured using CCD data and the 
most recent point is from our Robo-AO observation. The blue lines 
are the predicted change in separation (upper panel) and position 
angle (lower panel) assuming the star 8(^8 to the SW of KELT-8 is 
in the background with negligible proper motion and KELT-8 has 
a constant proper motion with the values given in Table The 
dashed black line is a model assuming that the two components are 
physically bound and the orbital motion of the pair is negligible 
during the span of observations. 


observations using the synthetic template technique. We 
artificially inject Gaussian noise into the spectra before 
sending them through the normal RV extraction pipeline. 
We use an exposure meter to expose to S/N=200 for 
all normal observations and the highest S/N ratio data- 
point shown in Figure has no artificial noise injected. 
The RV RMS monotonically decreases down to extremely 
low S/N ratio levels with a noise floor at just below 10 
ms“^ for high S/N ratio. Two other RV standard stars 
were also run through the synthetic template pipeline 
with S/N ratio artificially decreased to match that of 
our KELT-8 observations. In each case the synthetic 
template RV RMS was about 10 ms“^ (see Figure [^. 
We use this technique to extract the RVs from the ARF 
data for KELT-8 and achieve a median per-measurement 
precision of 11.8 ms“^. 

While our technique is unique for iodine-based RVs in 
the optical, several other very similar techniques have 
previously been de velop ed to extract high p recision RVs. 


Bailey et al. (2012) and Tanner et al. (2012) use interpo- 
lated model atmospheres tor their ternplates, but instead 
of an iodine cell they use telluric lines for their PSF and 
wavelength calibration. They were able to achieve ^50 
m/s precision in the near infrared with NIRSPEC on 
Keck. Other telluric-calibrated techniques have proven 
successful in the optical with HIRES to obtain velocities 


to 100 m/s precision (Chubak et al.|2012 ). Johnson etal. 
(20061 perturb an existing empirical template of a similar 
star m order to avoid the need to collect a unique tem¬ 
plate for each star observed. They were able to achieve 
a prescision of 3-5 ms“^ using these perturbed empirical 
templates. 



Fig. 7.— Synthetic template RV performance as a function of 
S/N ratio of the iodine observations. Artificial Gaussian noise 
was injected into the spectra of a well-known RV standard star 
(HD 9407) before running them through the RV extraction pipeline 
for all data points excluding the one at S/N=200. The synthetic 
template technique shows predictable and consistent performance 
down to extremely low S/N ratios and a model-limited noise floor 
of RilO ms“^. 


5. PLANETARY CHARACTERIZATION 
5.1. Global Fit 

We determine the system’s orbital and planetary prop¬ 
erties by performing a global analysis on the RVs, spec¬ 
troscopic stellar parameters, and follow mp light curves 
using a custom version of EXOEAST (Eastman et al. 
2013 1. Our global analysis technique is ve ry similar to 
that used for the previous KELT discoveries (jSiverd et al. 


2012 Beatty et al.|2012 [Pepper et al.|201^|Colfins et al. 

■ -■ 


2014 Bieryla et al.||2015|). The EXOFAST package uses 

a DE-MOMO fitting routing to extract the median pa¬ 
rameter values and 68% confidence intervals using all 
available data and these values are presented in Table 

a Constraints on the stellar parameters M* and i?* are 
so included using the Yonsei-Yale stellar models and 
spectroscopically derived values for and [Fe/H] as 
input priors. Before running the MCMC chains, a best- 
fit solution is found using an AMOEBA minimization 
routine and the photometric and RV uncertainties are 
scaled such that x^=l. The scaling factors for the pho¬ 
tometric data are listed in Table [l] and the scaling factor 
applied to the APE RV uncertainties is 1.44. The scaled 
uncertainties include both instrumental and astrophysi- 
cal noise sources and this ensures that the widths of fi¬ 
nal posterior properly represent the uncertainties in the 
model parameters. 

As described in Collins et al. (2014), we identify trend 
datasets that best improve the individual light curve fits 
using AIJ. These trend datasets are included as input 
to EXOEAST, along with the undetrended light curve 
data. As part of the global fit, EXOEAST simultane¬ 
ously finds the best transit model fit and removes linear 
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Fig. 8.— Synthetic template RV performance for three RV stan¬ 
dard stars. Artificial Gaussian noise was injected into all spectra 
for each star so that S/N=50 for each spectrum. This is approx¬ 
imately the same S/N ratio that was collected for each KELT-8 
observation. The RV RMS is about 10 ms“^ for each star. 

trends that are correlated with different combinations of 
airmass, time, FWHM, and CCD position. 

The light from the visual companion to KELT -8 was 
included in the photometric aperture of the GCorfini, 
and ZRO light curves. We corrected for this by measur¬ 
ing the flux ratio of the primary to secondary component 
in the well resolved MORC data. The flux ratios are 1.24 
and 1.44 in g' and i' band respectively. We then interpo¬ 
lated/extrapolated these flux ratios to the central wave¬ 
lengths of the V and Ic GCorfini and ZRO light curves 
and corrected the input light curves for the dilution be¬ 
fore fitting them in EXOFAST. We also ran the global 
fit excluding these two transits and found that all pa¬ 
rameters are consistent to well within 1 a. This gives us 
confidence that these blended light curves are not biasing 
the inferred parameters. 

We ran four permutations of the global fit using slightly 
different input priors and constraints. First, we as¬ 
sert constraints on the stellar mass and radius from the 
Yonsei-Yale isochrones and a uniform prior on log(/. We 
adopt the parameter values from this fit for all analy¬ 
sis and interpretation. Second, we performed the same 


fit using constraints from the Torres et al. (|2010 ) rela¬ 
tions instead. These two fits produced nearly identical 
results. Then we performed these same two global fits 
but imposed a prior on log 5 from the high-resolution 
spectroscopy. The logg derived from spectroscopy is 
slightly higher and only marginally consistent with the 
log 5 derived from transits. The results from all four fits 
are presented in Table The differences in parameter 
values from these four fits give an estimate of our sys¬ 
tematic uncertainty which appears to be similar to the 
statistical uncertainties for most parameters. This gives 
us confidence that our parameter uncertainties are prop¬ 
erly estimated and that our interpretations of physical 
parameters (such as the large planetary radius) are ro¬ 
bust. 

We measure an eccentricity for KELT- 8 b of e = 
0.035^0 025 which deviates from zero by only 1.4 a. The 
measurement of eccentricity is easily biased to artificially 
larger values and a significance of >2.5 a is generally ac- 
cepted as being require d to claim a non-circular orbit 
(Lucy & Sweeney 1971). Although we do not find any 
significant eccentricity in the system, moderate eccen¬ 
tricity is not unprecedented in other hot Jupiters (e.g. 
HAT-P-2b, XO-3b, WASP14b, HAT-P-31b |Bakos et al. 


2007 


Winn et ari|2009 Joshi et n)]|2009| [Kipping et al. 


20111 ! We allow for eccentricity to vary as a free pa- 


rameter in order to ensure that the errors on all other 
parameters are not underestimated. 

We searched for possible transit timing variations 
(TTVs) in the system by allowing the transit times for 
each of the follow-up light curves to vary. The ephemeris 
is constrained by the RV data and a prior imposed from 
the KELT discovery data. We find no evidence for sig¬ 
nificant transit timing variations. The ind ivid ual transit 
times are presented in Table and Figure [TOl 

After the completion of the global fit we then derive a 
linear ephemeris. The period is constrained by the RV 
data and comes from the global fit. We fit for To us¬ 
ing the transit times from the followup light curves. The 
best-fit linear ephemeris is listed in Table^ We find P = 
3.24406 ± 0.00016 days and Tq = 2456883.4803 ± 0.0007 
BJDtdb- The for the linear ephemeris fit is 7.4 
with 4 degrees of freedom which indicates that a lin¬ 
ear ephemeris adequately describes the observed transit 
times. 


5.2. False-Positive Assessment 

We calculate bisector spans (BS) for each of the APE 
spectra in order to investigate possible blend scenarios 
in which the measured RV variations are caused by line 
distortions that arise from an unresolved eclipsing binary 
instead of true reflex motion of the star. In the case 
that the measured Doppler shifts are caused by a second 
set of incompletely resolved stellar lines from a luminous 
companion we would expect to see BSs that vary in phase 
with the orbital signal and a correlati on in BS vs. RV 


(Santos et al. 

2002 

Torres et al. 2005). The amplitude 

would be simi 

lar to that of the velocities (Queloz et al. 

2001 

Mandushev et al.||2005|. 

Eo 

lowing the prescription of 'Lorres et al. (2005) and 

Torres et al. 

(2007 

), each API' spectrum was cross- 


correlated with a synthetic spectrum derived using the 
stellar properties determined from the high S/N ratio 
Keck spectrum (see Section]^. We analyzed the cross- 
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TABLE 4 

Median values and 68% confidence intervals for KELT-8b. 


Parameter 

Units 

Adopted Value 

Value 

Value 

Value 



(YY) 

(Torres) 

(YY±logp prior) 

(Torres±logp prior) 


Stellar Parameters: 

M* . Mass (Mq ). 

. Radius (Rq ). 

L*. Luminosity (Lq ). 

p*. Density (cgs). 

log p* .... Surface gravity (cgs). 

Teff. Effective temperature (K) 


-1 oil -1-0.078 

1 ooq”I”0.066 

i.Z99_Q 003 

1 1 oq+0.067 

t-to5_0 0g3 

1 91-1 +0.060 
i.Zii_o 059 

1 ny+O-l*! 

1 no+O-lS 

-l-ll^-0.12 

1 e:7n+0-099 
i.O/U_0 097 

1.577 ±0.096 

9 74+0-51 
^•^^-0.40 

2 79+0-54 

^•^^-0.32 

2 44+®-^‘l 
^•^^-0.31 

n *^<19+0-073 

U.9Dy_Q Qg7 

0 2n*:i+0 078 
U.9D9 _q 072 

0 4*30+0.076 
U.499_0 0g;^ 

U.49D_0 002 

4.078l°;®« 

4.U/O_0 0g;L 

4 121 +0-043 
^•1^1-0.041 

^.1ZO_0 042 

5754l®4 

5749 ± 59 

5748 ± 53 

5750 ± 59 


[Fe/H] ... 

Metallicity. 

0.272 ± 0.038 

0.270 ± 0.039 

0.270l°;®®« 

p. 9 - 7 .+0.040 
u.zm _0 039 

Planetary Parameters: 





e. 

Eccentricity. 

0.035lO-®®25 

0 037+0-057 

U.UO/ _o 027 

0.028lg;®o?? 

0 027+^ *^^^ 
u.uz/_o 020 

i-LJ*. 

Argument of periastron (degrees). 

85l|? 

or+82 

^^-93 

Q/-»+i 20 

®^-130 

or» + 120 

°^-130 

a . 

Semi-major axis (AU). 

0.0457i1S;°So84 

0.045991°;®®°?! 

U.U^O^Z _0 00032 

U.U^O/i _0 00075 

Mp . 

Mass (Mj). 

0.8671°;®®? 

0-8761^6? 

U.OOO _0 00 ;^ 

n SBQ+O-Osi 
U.8b9_0 000 

Rp . 

Radius ( Rj) . 

i. 8 b _0 ^0 

1 ftS+ 0-10 
J-.O »_0 j 0 

1.73 ±0.13 

1.74 ±0.13 

PP . 

Density (cgs). 

0 157+0-047 

O'-iO /_0 033 

0 153+0-049 

u.±oo_o 039 

u.zuo _0 039 

U.ZUO _0 040 

loggp ... 

Surface gravity. 

9 7'qo+ 0.072 
z. <yo _0 075 

9 ' 7 Q/:? + 0.076 
z./ou _0 080 

Z.OOU _0 003 

9 050 + 0.066 
Z.OOO_0_004 


Teq . Equilibrium temperature (K). 

©. Safronov number. 

{F) . Incident flux (10^ erg s“^ cm~^). 

RV Parameters: 

Tp . Time of periastron (BJDtdb-2450000) 

K . RV semi-amplitude (m/s). 

Mp /M+ .. Mass ratio. 

lAPF ■■■■ m/s. 

ecoscj*. 

esino;*. 

Primary Transit Parameters: 

Rp/R:^ .. Planet to star radius ratio. 

a/R* . Semi-major axis in stellar radii. 

i' . Inclination (degrees). 

b . Impact parameter. 

5. Transit depth. 

T. Ingress/egress duration (days). 

Ti 4 . Total duration (days). 

Uij . Linear Limb-darkening. 

U 21 . Quadratic Limb-darkening. 

Uisioang' Linear Limb-darkening. 

U 2 sioang' Quadratic Limb-darkening. 

'^iSioani' • Linear Limb-darkening. 

U 2 sioani' • Quadratic Limb-darkening. 

Uisioanr' Linear Limb-darkening. 

U 2 sioanr' Quadratic Limb-darkening. 

UiY . Linear Limb-darkening. 

U 2 V . Quadratic Limb-darkening. 

Linear Ephemeris from Follow-up Transits: 

P . Period (days). 

To. Time of transit (BJDxdb-2450000) ... 

Predicted Secondary Eclipse Parameters: 


1675^55 
0 n‘^'51 +0.0040 

U.U001_Q QQ37 

1 ' 70 -I- 0.27 
-*-• '°-0.22 


6870.47t°J® 

104.0 ± 6.4 

n 000^^^*^+0 000045 
(J.UUUD53_Q QQQQ 44 

6.1 ±4.9 

U.UUi_Q Q21 

0 012+0-003 

U.U1Z_Q 025 


16791®® 

U.UO^O_Q QQ3g 

1 RO+0-^0 

J-.oU_o 23 


6870.471®;™ 
104.0 ± 6.4 

U.UUUD/y_Q QQQQ43 

6.0 ±5.0 

0 qoi+ 0 025 
U.UUi_Q Q21 

0 014+0 071 
U.Ui4:_Q q27 


16291^^ 

O ( 3070 -1-0.0039 

U.U6<'5_q QQ35 
-*^•0^-0.18 


6869 . 0 II 2 

104.1 ± 6.5 

O OOORS7+0-000045 

U.0UUb57_Q QQQQ44 
6.7 ±4.9 

0 QQ1+0 022 

U.UU1_Q 

— 0 QQ1+0-029 

U.UUi_Q Q32 


16281 :^® 
o 0*^77+0-0039 

U.KJOI ^_0.0035 
1 59+0-10 


6869. 0 II 2 

104.4 ± 6.5 
0.000685 ± 0.000044 

0 001^®'®^^ 
U.UUi_g Qjg 

— 0 001^®'®^® 
U.UU 1 _Q Q 34 


0.1145 ±0.0026 

0.1146 ±0.0026 

0.1133 ±0.0025 

0.1132 ±0.0025 

9-yu_o 38 

5 07 + 0.39 

^•°'-0.42 


p* 90 + 0.35 
^•^'^-0.31 

82.651?;?! 

82.561?;?® 

83.36 ±0.69 

0*3 QO+0.69 
83.38_0 70 

Q 741 + 0.027 
'4* '^-L_o.o33 

0 743+0-027 

u./^o_o 033 

n 799+0.028 

u./zz_o 036 

n 799 + 0.029 

u./zz_o 037 

).01311± 0.00059 

0 0131 3+0-00061 
u.uioio _0 00059 

n m 90*3+0.00058 
U.Ui283_0 00057 

0.01282 ±0.0005'i 

0 Q 3 Q 3+0 0038 

U.U9U9_0 0035 

U.UOUO _0 0030 

O.O28U_0 0032 

n 0980+0-0034 

O.O28O_0 0032 

0 1444+®-^®^'^ 

u.i ^^^_0 0033 

0.1446 ±0.0034 

0 1423+®-®®^^ 
U.l^ZO _0 003;,^ 

0.1423 ±0.0032 

0 2961+®'^®^^ 
u.zyui _0 0033 

u.zyDO _0 0093 

n 9Q77+0.0089 
u.zy / / _o.oos 7 

0 9Q7K+0.0096 
u.zy/o _0 0092 

0 2793+*^-^*^^'^ 
u.z/yo _0 0049 

u.z/yu _0 0052 

n 97Q9+0.0042 

0.2/^82_0 0047 

n 97OQ+0.0046 
U.2/'83_0 005;^ 

0.607 ±0.015 

0 507+0-016 
U.OU /_0 0^5 

0.608 ±0.015 

0 . 6081 ?;??? 

n 1 Q^+o.oio 

U.184_0 0 ^^ 

n 1 RA+O-Oii 

O-lo4_0 0^2 

n 10 * 3 + 0.010 

O.i83_0 0 j^ 

0.1841?;??? 

n 01 70 + 0 .0097 
u.oi/o _0 0092 

n *^ 10 * 3 + 0.010 

u.oioo _0 0097 

0 3193+0-0094 

u.oiyo _0 009;,^ 

0 3192+0-010 
u.oiyz _0 0090 

n 97OQ+0.0046 
U.2/89_0 0052 

n 970^+0.0050 

0 . 2 /' 8 b _0 0050 

n 977Q+0.0045 
u.z/ <y_o.oo5i 

0 977Q+0.0049 
u.z/ <y_0.0055 

u.^iu_o on 

0.410 ±0.012 

0 411+0-012 
o.^ii_o.on 

0.411 ± 0.012 

n 9795+0.0062 
u.z/zo _0 0070 

n 9799+0.0066 

U.Z/ZZ _0 0074 

n 971 5+0.0062 

U.z/ io _0 0003 

0 971 7 + 0.0066 

U.Z/l^_o 0074 

0.4871?-®?? 

0.4881?;??? 

0.4881?;??? 

0.4881?;??? 

0-24461?;??™ 

0 2442 

u.z^^z _0 0093 

0.24351?;??™ 

0 9^*37+0.0082 

U.z^o <_0 0093 


3.24406 ± 0.00016 
6883.4803 ± 0.0007 


3.24406 ± 0.00016 
6883.4803 ± 0.0007 


3.24406 ± 0.00016 
6883.4804 ± 0.0007 


3.24406 ± 0.00016 
6883.4804 ± 0.0007 


Ts . 

Time of eclipse (BJDtdb-2450000) ... 

6868.8861?-??? 

6868.8871?-??? 

6868.8861®;°?? 

6872.1301?;??? 

bs . 

Impact parameter. 

0 753+0-097 

U. ^D9_0 057 

0-7691®;??! 

0.7181?;??? 

0.7161°;?®? 

'TS . 

Ingress/egress duration (days). 

U.U9Z/_0 0059 

U.U99Z_0 0005 

n 0977+0.0054 
V.VZi /_0.0042 

0 0275+0-0053 
U.U^/D_0 0042 

'7s'.14 • ■ • ■ 

Total duration (days). 

0 1435+0-0034 
U.149D_0 0049 

0 1431^+0-0036 
U.i4:C)O_0 0000 

0 1419+0-0033 

U.14iy_0 004;,^ 

n 1A1 8+0-0034 

U.1418_0 0042 
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Fig. 9.— Top: APF RVs and residuals for KELT-8. The best-fit 
model is shown in red. Bottom: RVs and bisector span measure¬ 
ments phase-folded to the best-fit linear ephemeris. The best-fit 
model is shown in red. The predicted Rossiter-McLaughlin effect 
assumes perfect spin-orbit alignment and it is not constrained by 
our data. 


TABLE 5 

Transit Times for KELT-8b. 


Epoch 

Tcf 

(BJDtdb) 

O-T 

(s) 

O-C 

(s) 

O-C 

(oTc) 

Telescope 

-4 

2456870.502251 

200 

-156.3 

-0.78 

ZRO 

-4 

2456870.509188 

203 

443.1 

2.18 

GCorfini 

-3 

2456873.746945 

114 

-101.5 

-0.89 

MORC 

-3 

2456873.749450 

137 

114.9 

0.83 

MORC 

5 

2456899.699249 

138 

-116.7 

-0.85 

WCO 

8 

2456909.433338 

177 

48.2 

0.27 

CROW 


15 


10 


'S 

u 0 

6 

-5 

-10 

Fig. 10.— Residuals of the transit times from the best-fit 
ephemeris for the followup light curves. The source of each transit 
time is denoted in the upper right. The dashed lines indicate the 
1 a uncertainty on the linear ephemeris. 

correlation function for each of the echelle orders between 
4260 and 5000 A in order to avoid the iodine lines, and 
telluric lines. By restricting the BS analysis to blue or¬ 


• ZRO 

4 GCorfini 
■ MORC 
A WCO 

* CROW 



-4 -2 0 2 4 6 8 

Epoch Number 


ders we reduce the effect of instrumental PSF variations 
caused by guiding errors because the seeing is degraded 
toward the blue which means that the slit is more evenly 
illuminated. We measure the asymmetry in the spectral 
line profile for each order by calculating the velocity at 
the midpoint of lines connecting the CCF at many dif¬ 
ferent fractional levels of CCF peak. The BS is then 
the difference in velocity between the 65th and 95th per¬ 
centile levels of the CCF. The bisector spans and error 
corresponding to a given observation are the mean and 
standard deviation on the mean over the 15 spectral or¬ 
ders analyzed. 

Our BS measuremen ts a nd errors are listed in Table 

and shown in Figure Figure P also shows the BS 
measurements as a function of orbitm phase. We find no 
statistically significant correlation of BS with RV. The 
Spearman rank correlation coefficient is -0.35 (p=0.24, 
Spearman ||1904|. 

'i'he fack of in-phase BS variations, and the fact that 
the stellar log 5 derived from spectroscopy and logg de¬ 
rived from the transit light curves are consistent to within 
1-CT lead us to conclude that the RV variations and transit 
signals are caused by a highly-inflated Jovian planet. 



Fig. 11.— BS measurements for the APF spectra used for radial 
velocity measurements. We find no significant correlation between 
RV and BS. 


6. DISCUSSION 

6 . 1 . Prospects for further characterization of the host 

star 

Further stellar characterization of KELT -8 is needed 
to better constrain the mass and radius of the star which 
will improve our knowledge of the planetary mass and ra¬ 
dius. Several more light curves with higher precision in 
which the two components are well-resolved would help 
to better constrain the stellar density. A more detailed 
analysis of higher S/N ratio spectra may help to nail 
down the spectroscopically derived parameters and allow 
for chem ical abundance measurements. GAIA (|de Brui-j 
jne ||2012 1 should provide a precise parallax to constrain 
the stella r luminosity and th us radius, and TESS and/or 
PLATO (Rauer et al. 2014 1 may provide sufficient pho¬ 
tometric precision to measure the stellar density using 
asteroseismology. 


6 . 2 . Comparative Planetology 
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KELT-8b has the 2nd largest radius of all known exo¬ 
planets, but has a mass slightly less than that of Jupiter. 
It lies well above the theo reti cal mass-radius curve for 
pure hydro gen (see Figure [T^. KELT-8 b is similar to 
WAS P-lb (ICollier Camero n et al. 20071, and HAT-P- 
41b (Hartrnan et al. 2012" It orbits a star more mas¬ 
sive and metal rich than the Sun (M*> 1 Mq and 

[Fe/H]> 0.1), the planet is highly irradiated (a < 0.05 
AU), and it is less massive than Jupiter but extremely 
inflated {Rp > 1.5 i?j). 

Next, we compare the stellar parameters of KELT- 
8 with other well-characterized hot Jupiter host stars. 
We use the exoplanets.org catalogue to select the stars 
that host planets with orbital periods shorter than 10 
days, and masses between 0.4 and 1.1 Mj. We plot 

these stars on an HR diagram in Figure We split 

this sample into two groups; one group containing the 
stars that host planets with radii greater than 1.5 i?j 
and the other group containing stars that host planets 
smaller than 1.5 i?j. KELT-8 is the coolest star to host 
a highly inflated hot Jupiter {Tes = 5754 K), but we find 
that all other stars that host similar highly inflated plan¬ 
ets tend to be hotter than wSOOO K. However, there are 
also plenty of stars hotter than 5900 K that host non- 
inflated hot Jupiters. We also compare stellar metalli¬ 
cally, V sin i, flux received by the planet, orbital eccentric¬ 
ity, stellar multiplicity, and spin-orbit misalignment, but 
find no clear differences between the hot stars that host 
“normal” hot Jupiters to those that host the highly in¬ 
flated hot Jupiters with radii > 1.5 i?j. However, we note 
that several of these host star parameters are poorly con¬ 
strained for many of the systems. Albedo would be an¬ 
other interesting parameter to compare for these planets 
but it has not been measured for the majority of planets. 
It is possible that the highly inflated planets have ab¬ 
normally low albedos and absorb more energy int o their 
atmosphere s at a given stellar insolation. TESS (Ricker 
et al.||2014 1 should provide light curves of sufficient pre¬ 
cision to measure secondary eclipses in reflected light for 
most or all of these systems. 


6.3. Mass-Radius predictions 

Several authors have identified empirical relations that 
attempt to predict the radius of a planet bas ed on a com¬ 
binati on of stellar and planetary properties. Enoch et al. 
(2012) established a relation that predicts planetary ra¬ 
dius based only on the equilibrium temperature and or¬ 
bital semi-major axis. KELT-8 falls i nto the Jupiter mass 
bin as defined by Enoch et al. (20121 and the relation for 
that mass regime predicts a radius of 1.43 Rj which is 
far too small relative to our measured valu e. The most 
recent empirical relation was presented by [Weiss et al.| 
(2013), which uses the incid ent flux and planet ary mass 
as dependent variables. The Weiss et al. (2013) relation 
predicts a radius of 1.63 R.t f or KELl'-8. I t one takes 
into account the scatter in the Weiss et al. (2013) rela¬ 
tion (0.1 R]), and the error on our radius measurement 
the predicted and measured values agree to within 1.2 ct. 


6.4. Prospects for atmospheric characterziation 

KELT-8 joins a rare breed of highly inflated hot 
Jupiters orbiting relatively bright stars. The extended 
atmosphere, combined with the large transit depth and 



Fig. 12.— Mass radius diagram of all confirmed transiting ex- 
oplanet^ (red circles). KELT-8b is the annotated black square. 
Solar system planets are represented by green triangles. Model 
mass-radius relationships for idealized planets consisting of pure 
hydrogen, water, ro ck (Mg 2 Si 04 ), and iron are shown as blue lines 
ijSeager et al.|2007|). Only WASP-17b has a radius larger than the 
best-nt radius lor KELT-8b. 

‘‘based on a 2/17/2015 query of exoplanets.org 


bright apparent magnitude makes KELT-8 one of the 
best targets for transmission spectroscopy. In addi¬ 
tion, the large planetary radius and high equilibrium 
temperature give rise to a significant secondary eclipse 
depth. We calculate an expected eclipse depth of «1.3 
mmag at 3.6 microns by taking the ratio of the black- 
body emission from the planet to that of the star and 
multiplying by A secondary eclipse of this 

depth would be easily accessible by Spitzer. Combining 
Teq = 1675 from Tablej^with the planetary surface grav¬ 
ity of \og gp = 2.793, we find that the scale height of a 
H2-dominated atmosphere {nm = 2) on KELT-8 would 
he H = 1113 km. This large scale height implies that the 
amplitude of a trans mission spectr oscopy signal would be 
^^£ 2 ^ = 0.22 mmag ( Winn||2010 assuming H « Rp). 
All of these calculations assume perfect heat redistribu¬ 
tion and zero albedo for the planet. Studies of the at¬ 
mospheric composition or temperature profile of KELT-8 
may provide an explanation for its highly inflated radius. 


6.5. Irradiation History 

We have shown above that KELT-8b is a highly in¬ 
flated planet, joining the ranks of other hot Jupiters that 
manifest radii much larger than predicted by standard 
models for non-irrad iated objects with Jovia n masses. 
Several authors (e.g., Demory & Seager 2011) have sug¬ 
gested an empirical insolation threshold (« 2 x 10® erg 
s“^ cm“^) above which hot Jupiters exhibit increasing 
amounts of radius inflation. KELT-8b clearly lies above 
this threshold, with a current estimated insolation of 
1.78lo;22 10® erg cm“®, and therefore its currently 
large inflated radius is not surprising. At the same time, 
the KELT-8 host star is found to currently be in a very 
rapid state of evolution, such that its radius is rapidly ex¬ 
panding as the star crosses the Hertzsprung gap toward 
the red giant branch. This means that the star’s surface 
is rapidly encroaching on the planet, which presumably 
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Fig. 13.— HR diagram showing well-characterized host stars of 
hot Jupiters. The red squares indicate stars that host planets with 
radii greater than 1.5 Rj. The green star marks the position of 
KELT-8, and the black circles are all other systems. 


is rapidly driving up the planet’s insolation and also the 
rate of any tidal interactions between the planet and the 
star. 

Therefore it is interesting to consider two questions. 
First, has KELT-8b’s incident radiation from its host 
star been below the empirical radius inflation thresh¬ 
old in the past? If KELT-8b’s insolation only recently 
exceeded the inflation threshold, the system could then 
serve as an empirical testbed for the different timescales 


predicted by different inflation mechanisins (see , e.g., As- 
sef et al.||2009t iSpiegel & Madhusudhan||2012|). Second, 
what is the expected fate of the KEL'i'-8b planet given 
the increasingly strong tidal interactions it is experienc¬ 
ing with its encroaching host star? _ 

To i nvestigate these questions, we follow [Penev et ah] 
(2014) to simulate the reverse and forward evolution of 
the star-planet swtem, using the measured parameters 
listed in Tables ^ and as the present-day boundary 
conditions. This analysis is not intended to examine 
any type of planet-planet or planet-disk migration ef¬ 
fects. Rather, it is a way to investigate (1) the change in 
insolation of the planet over time due to the changing lu¬ 
minosity of the star and changing star-planet separation, 
and (2) the change in the planet’s orbital semi-major 
axis due to the changing tidal torque as the star-planet 
separation changes with the evolving stellar radius. We 
include the evolution of the star, assumed to follow the 
Yonsei-Yale stellar model with mass and metallicity as 
in Table For simplicity we assume that the stellar 
rotation is negligible and treat the star as a solid body. 
We also assume a circular orbit aligned with the stellar 
equator throughout the full analysis. The results of our 
simulations are shown in Figure We tested a range 
of values for the tidal quality factor of the star from 
log <5* = 5 to log (5* = 7 (assuming a constant phase lag 
between the tidal bulge and the star-planet direction). 
Q* is defined as the tidal quality factor divided by the 
Love number = Q*/k 2 ). We find that although for 
certain values of Q'^ the planet has moved substantially 
closer to its host during the past Gyr, in all cases the 
planet has always received more than enough flux from 
its host to keep the planet ir radiated beyond the inso - 
lation threshold identified by Demory & Seager (2011), 
except perhaps during the pre-main-sequence (prior to 



0.1 1.0 
Age [Gyr] 



Fig. 14.— Top: Irradiation history of KELT-8b. The insola¬ 
tion received by the planet is well above th e empirical inflation 
irradiation threshold (|Demory Sz Seager|2011|l for the entire main- 
sequence existence of ttie star except m ttie case of log Q'^ = 5 in 
the very early stages of stellar evolution. Bottom: Orbital semi¬ 
major axis history of KELT-8b. The planet’s semi-major axis is 
rapidly decreasing as the star evolves off the main sequence. It 
appears unlikely that KELT-8b will survive past the star’s current 
subgiant phase. 

an age of ~100 Myr). 

Interestingly, the currently rapid evolution of the star 
suggests a concomitant rapid in-spiral of the planet over 
the next few 100 Myr, unless the stellar Q'^ is large. This 
planet therefore does not appear destined to survive be¬ 
yond the star’s current subgiant phase. As additional 
systems like KELT-8b are discovered and their evolution 
investigated in detail, it will be interesting to examine 
the statistics of planet survival and to compare the se 
to predictions such as those shown here in Figure [14] to 
constrain mechanisms of planet-star interaction generally 
and the values of Q'^ specifically. 

7. SUMMARY 

We announce the discovery of the highly inflated hot 
Jupiter, KELT-8b. This planet was initially discovered 
in KELT photometry, then confirmed via high-precision 
followup light curves and RVs. We also present adaptive 
optics imaging of KELT-8 from the Robo-AO system. 
Astrometry measurements over the past 110 years com¬ 
bined with high resolution spectra of both components 
firmly establishes the visual companion at 8?8 separation 
as an unrelated background star. 
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We develop a new technique to extract reliable RVs 
from noisy data that saves a significant amount of tele¬ 
scope time. In the case of RV follow-up for the confirma¬ 
tion or mass measurements of transiting planets where 
the ephemeris is known and only a few well-timed RV 
measurements are needed this technique can save 30-50% 
of the time by avoiding the need to collect a high quality 
iodine-free template observation. In addition, lower S/N 
ratio is required for the RV measurement observations 
(S/N«50) through the iodine cell because one compo¬ 
nent of our model is now completely noise free. Normally, 
we collect S/N=200 for our precision RV exposures using 
the standard observed template technique. 

KELT-8b has one of the largest radii of any known 
transiting planet with Rp = I.86tQ;jg i?j. It joins a 
small and interesting class of highly inflated hot Jupiters 
orbiting stars slightly more massive than the Sun. 
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